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One main goal of quantum control is to steer a quantum system toward an expected 
state or dynamics. For measurement-induced quantum control, measurements serve as the 
only control, which is like the cases in quantum Zeno and anti-zeno effects. In this paper, 
this scenario will be investigated in a general A^-level quantum system. It is proved that, 
when the intial and expected states are both pure, the control space could be reduced to 
an effective subspace spanned by these two states only. This result will greatly simplify the 
measurement-induced control strategy of an TV-level quantum system. 
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I. INTRODUCTION 

Quantum control has drawn much attention both theo¬ 
retically and experimentally[ll-[^ due to its potential ap¬ 
plications. The goal of quantum control is usually to 
steer a quantum system into an expected quantum state 
via external fieldsHowever, the control means are 
not limited to introducing control fields. With rapid de¬ 
velopment of the research field, various approaches were 
proposed to play the role of control^. The quantum 
measurement, as one of the necessary criteria for a quan¬ 
tum computer, does not only serve as the read-out but 
also could be employed as the control^. 

The idea of utilizing measurements to drive quantum 
systems origins from the famous quantum Zeno effect 
which claims that repeated frequent measurement of an 
unstable state will keep it unchangedjl, [^. The coun¬ 
terpart of Zeno effect, namely anti-Zeno effect , was 
also discussed [13, [HI , and both of them were observed 
in an unstable systemfl^. In contrast to control field 
schemes, measurement-based quantum control will in¬ 
troduce irreversible decoherence inevitably. By combin¬ 
ing with coherent control, it can be applied to closed- 
loop molecular states control [T^. quantum dynamics 
engineering[l^, dephasin g d ecoherence controlflsl and 
remote state preparation 1 1 6| . Controllability analysis 
of quantum systems shows that measurements could 
transform a non complelely controllable system into a 
controllable one while unitary evolution only could not 
achieve this task[T3, Si- System controllability un¬ 
der Kraus-map, which includes non-projective measure¬ 
ments, was also studied [l3,113 ■ More extreme cases with 
measurement-only control schemes have been studied. In 
the absence of active coherent controls, projective mea¬ 
surements were introduced to map an unknown mixed 
state onto a known expected pure state while only two 
noncommuting observables are available!^. Optimized 
measurement se quen ce was obtained for two-level quan¬ 
tum systems!^ The continuous measurement ex¬ 
tension was investigated and it shows that quantum anti- 
Zeno effect could be recovered in the limit of continuous 


cases [3. [ 23 ]^ . Schemes for state engineering using mea¬ 
surements and fixed dynamics were also peoposed(T7|. 
Recently, a novel quantum state preparing protocol was 
presented in which only a restricted set of measurements 
is reauiredf^. In this protocol, previous measurement 
outcome is feedback as the control, and it shows that ar¬ 
bitrary state could indeed be prepared. This feedback 
idea was also introduced to state manipulation!^. 


These above works investigated both two-leve l iTa. 1^ . 
US, and Y-level systems [1, [3, [13, IS HH I27l.l3r)l|. 

In this paper, we will investigate how to control TV-level 
systems with arbitrary projective measurements only, 
which is a direct extension of our previous works [SOj- 
However, powerful tools in two-level context, i.e. Pauli 
matrics and Stokes vector, are no longer valid in TV-level 
systems with TV > 2. In Ref(^. a variational method 
over unitary group was introduced to obtain the neces¬ 
sary condition for optimal measurement sequence. This 
condition was carefully analyzed for pure initial and ex¬ 
pected states with a single measurement. This will be 
extended to arbitrary finite times of projective measure¬ 
ments for pure state-to-state control. We find that the 
conrol space will be reduced to an effective subspace 
spanned by the initial and expected states only. The 
transition cases through intermedia levels are excluded 
and the two-level modelwill be recovered. 


The rest of the paper is organized as follows. In Sec. 
H, we introduce some basic notations, the variational 
method over unitary group and the necessary condition 
for an optimal measurement sequence under non-selective 
assumption. In Sec. HI, we investigate an orthogonal 
state-to-state control case, and derive the effective sub¬ 
space for optimal control. These cases will be generalized 
to ones of non-orthognonal pure states in Sec. IV. Dis¬ 
cussion and conclusion are presented in Sec.V. 
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II. NECESSARY CONDITION FOR OPTIMAL 
STATE TRANSFER 

The measurement-based optimal control case for a two- 
level quantum system was solved by using Pauli matrices, 
and the optimal projectors were also obtainedHJ. Since 
this approach seems invalid for general TV-level systems, 
we introduce a virational method and derive a necessary 
condition for the optimal measurement sequence [13 • In 
this section, some basic concepts and notations are firstly 
be introduced. Then, the virational method and the nec¬ 
essary condition for the optimal measurement sequence, 
which appears as a chain of equalities, are presented . 
These tools will be utilized to studyA-level cases in the 
next section. 

Consideing an iV-level quantum system whose states is 
characterized by a density matrix p, the state after a non- 
selective measurement A4 becomes A4 (p) = PipPi 

where {Pi} is a set of rank-1 projectors satisfying the 
equations PiPj = SijPi and Pi = I. Because of these 
properties, we could always find a suitable unitary matrix 
U to diagonalize these projectors, i.e. Pi = U\i){i\W 
for i = I,-- - ,N, where |f) indicates a specific level of 
the system. With these notations, the non-selective mea¬ 
surement process could be recast as 

N 

M{U)p = '£u\z){i\U^pU\i)m^ ( 1 ) 

i=l 

which contains three successive operations of p: (1) ro¬ 
tate the basis (a unitrary transfermation) (2) set non¬ 
diagonal elements to be zero (3) rotate the basis inversely. 

Limited by the measurement techniques nowadays, it 
is difficult to performe too many projective measure¬ 
ments within finite time interval. Assuming that at most 
m-times measurement could be employed in this inter¬ 
val, the final state p™ after these measurements reads 
p"* = M{Um) ■ • ■ M{Ui)p. We neglect free evolution of 
the system because this could be included via a picture 
transformation [^. The control objective J, which 
quantihes this control effect, is defined as the overlap be¬ 
tween p™ and some expected final state O'. 

J = Tr{p^0)=Tr{{MiUm)---M{Ui)p)0). (2) 

The necessary condition for optimal control sequence 
is obtained by considering the first derivative of J with 
respect to variables {Ui}. A variational analysis method 
is introduced to parametrize the neighborhood of a uni¬ 
tary transformation and the necessary condition presents 
as a chain of equalities: 


[p,M{U^)...M{U^)0] 


=[M[Uk)...M{Ui)p,M{Uk+i)...M{Um)e] 


=[M{Um)-M{U^)p,0]. 


( 3 ) 


Detailed mathematical derivation of the above equalities 
and some special cases for single measurement was shown 
in a, but general cases with m-times measurements are 
still intractable. The rest of this paper will focus on this 
general scenario. 


III. EFFECTIVE SUBSPACE FOR OPTIMAL 
ORTHOGONAL PURE STATES TRANSFER 

In this section, the effective subspace for orthogonal 
pure state-to-state optimal measurement control will be 
investigated. Assuming that the initial and expected 
states are p = |1) (1| and 0 = |2) (2|, respectively, and all 
states are represented in the basis {|i),i = 1,2, ....N}. 
Due to the special form of p, [p, M.(Ui)...M{Um)0] be¬ 
comes a matrix with non-zero elements only in the 
first row and first column . By the same reason, 
\M{Um)---M.{Ui)p, 0] becomes a matrix with non-zero el¬ 
ements only in the second row and second column. Com¬ 
bining these two forms, we have 


[p,M[U^)...M{Um)e] 


= [M[Um)-M{Ui)p,0] 


o 

where c is a complex number and O is a zero matrix of 
suitable size. The case of c = 0 is dropped henceforth, 
since it corresponds to a trivial control which can not be 
the optimal one. 

To clarify the discussion below, some notations are in¬ 
troduced here: /" denotes n x n identity matrix; U{n) 
denotes n-dimensional unitary group and f/” represents 
some elements in U(n)', O" is n x n zero matrix; Din) is 
the set of n X n density matrices, i.e. positive semidefinite 
Hermitian matrices. With these notations, we will first 
investigate the symmetry property of {Ui} and propose 
lemma. 1. 

Lemma. 1 \i {Ui,U 2 , ■ ■ ■, Um} is an optimal control se¬ 
quence, then {UsUi,UsU 2 , ■ ■ ■, UgUm} could also achieve 
this optimum, where Ug = © Uarb- Here, "01 , 

02 are arbitrary phases and Uarb S U{N — 2). 

The proof of lemma. 1 is straightforward. Assuming 
that {U{Ui,UgU 2 , ■ ■ ■ ,U}{^Um} could also achieve this 
optimum, then we have sufficient conditions for optimal 
control: {i)Ul^Ul'^^ = I ioi I = 1, ... ,m — 1 and (ii) p 
and 0 are invariant under the action of U{ and C/™, re¬ 
spectively. Consequently, we find that U{ = ■ ■ ■ = U^ 
and U{ = U^ = © Uarb- 

According to lemma. 1, we could always find proper 0i 
or 02 to eliminate the phase in c and make it real. Thus, 
c will be treated as a real number hereafter. 

Now, we investigate the first two parts of Eq. O and 


0 c 
-c* 0 
O 
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get 


[p,M{Ui)t] = C (5) 

[M{Ui)p,t] = C (6) 


in which r is defined as Multipling 

Ul and Ui from left and right in Eqs.® and m, and 
taking the {p, q) matrix element of both sides , we have 


(UItUi)pp^ 

= {iUl)lpiU^) 2 q - iU,r^piU^),q) C (7) 

= ((t/i)tp(C^i) 2 , - (t/i);p(C/i)i,) c (8) 


Obviously, Eq.Q implies that if (C/i)ip ^ 0 and 
{Ui)iq ^ 0 then (C/|rf/i)gg — {uItUi)pp could be 

expressed by {(t/i)ii} and {(t/i) 2 i} explicitly. Fur¬ 
ther, we can write the objective J as Tr{M.{Ui)pT) = 
Yjt{Ui)ii{Ui)li{UlTUi)ii. Note that {{UlTUi)ii} could 
be reconstructed by {(C/i)ii} and {(17i)2i} up to a con¬ 
stant, and if we fix this constant, which corresponds 
to solving Eq.([7]) and Eq.® in a special case, then 
(C/|rf/i)ii will only depend on {(C7i)ii} and {{Ui) 2 i}, and 
so does J. This observation will give us the explicit form 
of r. 

Using lemma. 1, we could diagonize the submatrix 
spanned by the last N — 2 rows and columns in r, and 
then investigate Since (C/|t17i)ii only de¬ 
pends on {{Ui)ii} and {{Ui) 2 i}, we have = 


3H,J()J 1 Again, by lemma. 1, we can as¬ 

sign some Uarb to make {(C/i)ii} real for i = 3,... ,N. 
This will simplify our calculation below. Remember that 
there is a normalization constraint, i.e. = 

1, which implies that {(t/i)ii} are not independent in 
(C//rl7i)ii. Due this constraint, we treat (C/i)Ari as a 
non-independent variable and have 


d{Ui)3i 

_ d {{Ui)iiT 3 i + {Ui)21T32 + C.c.) (Uljai 

9(Ui)31 

d {{Ui)iiTni + {Ui) 21TN2 + c.c.) {Ui)n1 

d{Ui)3i 

^ ((^l)3l(^l)31'r33 + (Ul)^i(Ul)jVl'rAfjv) 
d{Ui)3i 

=0 (9) 


Since (C/i)ii, (Ui)ii, (t/i) 2 i, (^ 1)21 are indepentdent, 
three parts in Eq.® gives T 13 = T 31 = r 23 = T 32 = 0 , 
tni = tm 2 = Tin = T 2 n = 0 and T 33 = tnn- By the 
same reason, we also have tu = T 2 i = 0 and th = tnn if 

we investigate for i = 3,..., N. This implies 

r £ D(2) © where d is a normalized factor . 


Next, we start from this special form of r to explore 
what kind of {Ui} could achieve the optimum. We define 
T = Ts © dl^~‘^ where Ts £ D{2) and write M.{Ui)p 


in a block form, i.e. M.{Ui)p = 


in which 


^11 ^12 
A21 ^22 

All and ^22 are 2x2 and {N — 2) x {N — 2) matrices 
respectively. According to Eq.®, we get {dl"^ —Ts)Ai 2 = 
O whose solutions have two different branches: (a) Ts 
has no eigenvalue of d or (b) Ts has an eigenvalue of 
d, and these two branches indicate A 12 is zero or not, 
respectively. 

For case (a), we have ^12 = O and M.{Ui)p = 


^11 O 
O A 22 


To acquire this form of A4{Ui)p, Ui must 


be (a.l) Ul £ {U{2)®U{N -2)) x (1 © C/(A^ - 1)) or 
(a.2) Ul £ U{2)Q)U{N — 2). Note that these two different 
cases origin from whether An and ^22 have degenerated 
eigenvalue. 

For (a.l), we define Ui = ViWi in which Vi £ U(2) (B 
U{N — 2) and lUi £ 1 © 17(7V — 1). Since the rotation Vi 
makes p and t belong to D{ 2 )(B 0 ^~‘^ and D{ 2 )(BdI^~^ 
respectively, a further rotation Wi makes some diagonal 
elements, except the first one, in Vi pVi equal, which are 
the degenerated eigenvalue in An and ^ 22 - Wi will also 
make the diagonal elements in the correaspnding sub¬ 
space of VitVi equal. This is due to the fact that this 
two subspaces, i.e. diag{a, 0,..., 0) and diag{b, d,... ,d) 
are essentially of the same form of adiag{a, 0,..., 0) + 
Pdiag(l, !,...,!) which will have same diagonal elements 
under the act of Wi. This property also implies that 
Wi must act on the whole (TV — l)-dimensional subspace 
to achieve the maximum. Define the rotation with this 
property by Ue{N), we get Ui £ (17(2) © 17(TV — 2)) x 
(1 © 17e(TV — 1)). Now, we have U 2 £ U{2) © 

[/(TV - 2) from (a) and Ui £ ([/(2) © [/(TV - 2)) x 
(1 © Ue{N - 1)). This implies M{U 2 )M{Ui)p £ D{2) © 
j^tjN -2 j^Q degenerated eigenvalue exists be¬ 

tween these two subspaces. According to Eq. (|T|), 
it indicates that M{U 3 )...M{Um)d £ D{2) © dl^~^. 
Again, we have U 3 belongs to U{2) © [T(TV — 2) or 
([/(2) © U{N - 2)) X (1 © Ue{N - 1)). The latter one, 
i.e. U 3 £ ([T(2) © [/(TV - 2)) x (1 © Ue{N - 1)), implies 
that two eigenvalues in A 4 {U 2 )...Ai{Um) 0 , which are not 
d, are both larger or smaller than d. If they are both 
larger than d, then Ui will belong to U{2) © [/(TV — 2) 
but not ([/(2) © [/(TV — 2)) X (1 © [/g (TV — 1)). If they are 
both smaller than d, then the objective could not exceed 
d which is less than 1/2, and this is contrary to the known 
optimumjH, [ 2 ^. Thus, we have U 3 £ [/(2) © [/(TV — 2). 
This procedure could be operated in an iterative fashion, 
and we finally have Ui £ U (2)©[/(TV—2) for T = 2,..., TV. 
Note that if we choose Ui = Vi and discard Wi, then J 
will be larger than the one when we have Ui = Vi x Wi. 
Thus we have [Ti £ [/(2) © [/(TV — 2) for T = 1,..., TV. 

If we have (a.2), then both Ui and U 2 belong to U (2) © 
[/(TV — 2), and Eq.® implies that M{U 3 )...M{Um)d £ 
D{2) © D{N — 2). By considering the special form 
of T, we also have M.{U 3 )...M.{Um)(^ £ D{2) © dl^~'^ 
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which indicates that U 3 £ C/(2) 0 U{N — 2) or U 3 £ 
(C/(2) © t/(7V — 2)) X (1 © t/(7V — 1)). For the former 
one, we could consequently investigate C /4 which revives 
the discussion of (a. 2 ). And the later one is just case(b) 
when we treat C /3 in (a.2) as C /2 in case (b). We will also 
show that Ui in case(b) also belongs to C/(2) © U{N — 2) 
as C /2 in (a.2). So, this situasion, i.e. (a.2), is equivalent 
to case(b). 

For case (b), the optimal control between r and p 
is equivalent to the optimal control between two non- 
orthogonal pure states. This could be seen when we di¬ 
agonalize T and g^ diag[a, d,... ,d]. The solution of this 
case is obtained!^ and we have Ui £ C/(2) © U{N — 2). 
Define U 2 = V 2 W 2 where V 2 £ U{2) ® U{N — 2), 
1^2 £ 1 © U^{N - 1 ), and 

X = V^M{U3)...M{Um)0V2 

rp ^ 21 ,up 

-Dll D 

— ^ 21 ,down 

_ ^ 21 .up ^ 21 ,down ^22 

where B 21 .UP -^ 21 ,rfou/n ^'I’e {N — 2) dimensional 

row vectors. Since arbitrary U G P ® U{N — 2) 

will not change the diagonal elements in V 2 M.{Ui)pV 2 
and a successive IF 2 makes diagonal elements, except 
the first one, in V 2 M(Ui)pV 2 equal, we find that 
B 22 = dl^~‘^ and B 2 i,down = 0. Define k = 
ZZi li) (il W^V^M{Ui)pV 2 W 2 \i) {i\ and Eq.® 
gives Kx ~ = F 2 ^C'V 2 which implies i/ 21 ,up = 0. Thus 

we have M{U 3 )...M{Ura)d £ D{2) © d"I^~^ and again, 
C /3 £ C/(2) © U{N - 2) or C /3 £ (C/(2) © U{N - 2)) x 
(1 © C/e(Af — 1)). The latter one should be discarded 
since it will not achieve the maximum when C /2 £ 
(C/(2) © U{N — 2)) X (1 © U{N — 1)), and finally we have 
C /3 £ C/(2) © U{N — 2). Following the analysis of (a.l) 
above, we also have Ui £ C/(2)©C/ {N—2) ioi i = 1,..., N. 

To conclude, for orthogonal pure states control induced 
by projective measurements only, the control space will 
be reduced to a two-dimensional effective subspace which 
spanned by the initial and expected states only. 


IV. GENERALIZATION TO 
NON-ORTHOGONAL PURE STATES 


The analysis in Sec.Ill could be generalized to cases 
with non-orthogonal pure states, and Eq.® is still valid. 
Define these two pure states by p = |1) (1| and 9 = 
(a |1) + /? |2)) {a* (1| + (3* (2|). Since p = |1) (1|, we have 

fo c si 


[p,M(C/i)...AI(C/™)0] 


c* 0 O 
st O O 


where c is a com¬ 


plex number and s is a {N — 2)-dimensional row vector. 
On the other hand, if we asuume that 9 could be diago¬ 
nalized by M, i.e. u9v3 = |1) (1|, then 


(m©/^“ 2 ) 


0 c s 
c* 0 O 
s+ O O 




= {u® I^-^)[M{Um)-M{Ui)p,9]{u^ ® I^-^) 

((«©/^-Xnt©j^-2)),|i)(i|] (10) 


The right hand of Ea. dTOll indicates that it will be a 
matrix with non-zero elements in the first row and col- 

0 c s 


umn only, thus {u © 


c* 0 O 
st O O 


(ui- © /^- 2 ) = 


'0 c' s'" 
c'* 0 O 
s't O O 


Since u is nontrivial, we must have s = O 


and Eq.® is recovered. The rest of the proof for non- 
orthogonal pure states is the same as that in Sec.III. We 
could also conclude that the effective subspace will be 
reduced to the one spanned by |I) and | 2 ). 


V. DISCUSSION AND CONCLUSION 


In summary, we discuss the use of projective measure¬ 
ments to control an V-level quantum system. Previous 
works of two-level scenario are generalized assuming that 
the measurement time is limited. It is proved that for 
pure initial and expected states, the control space will be 
reduced to an effective subspace spanned by these two 
states. The intuition for this result is straightforward, 
no extra level is needed for optimal control. The situ¬ 
ation for mixed states is much complicated. Numerical 
calculations show that this conclusion still works if “pure 
states” are replaced by “mixed states”, and its rigorous 
proof will be further investigated. 


ACKNOWLEDGMENTS 


This work is supported by the National Natural Sci¬ 
ence Foundation of China under Grant Nos. 61403362, 
61473199, 61203061, 61374091 and 11574295. 


[1] H. Rabitz, Science 288, 824 (2000), ISSN 1095-9203, [2] M. Shapiro and P. Brumer, Principles of the quantum 

URL http://dx.doi.org/10.1126/science.288.5467. control of molecular processes, vol. 1 (2003). 

824. 













5 


[3] D. d’Alessandro, Introduction to quantum control and dy¬ 
namics (CRC press, 2007). 

[4] 1. Petersen and D. Dong, lET Control Theory & Appli¬ 
cations 4, 2651 (2010), ISSN 1751-8652, URL http:// 
dx.doi.org/10.1049/let-eta.2009.0508. 

[5] A. Pechen, C. Brif, R. Wu, R. Chakrabarti, and H. Ra¬ 
bitz, Phys. Rev. A 82 (2010), ISSN 1094-1622, URL 
http://dx.doi.org/10.1103/PhysRevA.82.030101. 

[6] G. M. Huang, Journal of Mathematical Physics 24, 2608 
(1983), ISSN 0022-2488, URL http://dx.doi.org/10. 
1063/1.525634. 

[7] H. M. Wiseman and G. J. Milburn, Quantum measure¬ 
ment and eontrol (Cambridge University Press, 2009). 

[8] B. Misra and E. C. G. Sudarshan, Journal of Mathe¬ 
matical Physics 18, 756 (1977), ISSN 0022-2488, URL 
http://dx.doi.org/10.1063/1.523304. 

[9] W. Itano, D. Heinzen, J. Bollinger, and D. Wineland, 
Phys. Rev. A 41, 2295 (1990), ISSN 1094-1622, URL 
http://dx.doi.org/10.1103/PhysRevA.41.2295. 

[10] A. Balachandran and S. Roy, Physical Review Letters 
84, 4019 (2000), ISSN 1079-7114, URL http://dx.doi. 
org/10.1103/PhysRevLett.84.4019. 

[11] M. Lewenstein and K. Rz^zewski, Phys. Rev. A 61 
(2000), ISSN 1094-1622, URL http://dx.doi.org/10. 
1103/PhysRevA.61.022105. 

[12] M. C. Fischer, B. Gutierrez-Medina, and M. G. Raizen, 
Physical Review Letters 87 (2001), ISSN 1079-7114, URL 
http://dx.doi.org/10.1103/PhysRevLett.87.040402. 

[13] R. Judson and H. Rabitz, Physical Review Letters 68, 
1500 (1992), ISSN 0031-9007, URL http://dx.doi.org/ 
10.1103/PhysRevLett.68.1500. 

[14] S. Lloyd and L. Viola, Phys. Rev. A 65 (2001), 
ISSN 1094-1622, URL http://dx.doi.org/10.1103/ 
PhysRevA.65.010101. 

[15] M. Zhang, H. Dai, Z. Xi, H. Xie, and D. Hu, Phys. Rev. 
A 76 (2007), ISSN 1094-1622, URL http://dx.doi.org/ 
10.1103/PhysRevA.76.042335. 

[16] A. Mandilara and J. Clark, Phys. Rev. A 71 (2005), 
ISSN 1094-1622, URL http://dx.doi.org/10.1103/ 
PhysRevA.71.013406. 

[17] R. Vilela Mendes and V. 1. Man’ko, Phys. Rev. A 67 
(2003), ISSN 1094-1622, URL http://dx.doi.org/10. 
1103/PhysRevA.67.053404. 

[18] D. Sugny and C. Kontz, Phys. Rev. A 77 (2008), 
ISSN 1094-1622, URL http://dx.doi.org/10.1103/ 


PhysRevA.77.063420. 

[19] R. Wu, A. Pechen, C. Brif, and H. Rabitz, J. Phys. A: 
Math. Theor. 40, 5681 (2007), ISSN 1751-8121, URL 
http://dx.doi.org/10.1088/1751-8113/40/21/015. 

[20] A. Pechen, D. Prokhorenko, R. Wu, and H. Rab¬ 
itz, J. Phys. A: Math. Theor. 41, 045205 (2008), 
ISSN 1751-8121, URL http://dx.doi.org/10.1088/ 
1751-8113/41/4/045205. 

[21] L. Roa, A. Delgado, M. L. Ladron de Guevara, and A. B. 
Klimov, Phys. Rev. A 73 (2006), ISSN 1094-1622, URL 
http://dx.doi.org/10.1103/PhysRevA.73.012322. 

[22] A. Pechen, N. Il’in, F. Shuang, and H. Rabitz, Phys. Rev. 
A 74 (2006), ISSN 1094-1622, URL http: //dx.doi.org/ 
10.1103/PhysRevA.74.052102. 

[23] F. Shuang, M. Zhou, A. Pechen, R. Wu, O. Shir, and 
H. Rabitz, Phys. Rev. A 78 (2008), ISSN 1094-1622, URL 
http://dx.doi.org/10.1103/PhysRevA.78.063422. 

[24] F. Shuang, A. Pechen, T.-S. Ho, and H. Rabitz, 
The Journal of Chemical Physics 126, 134303 (2007), 
ISSN 0021-9606, URL http://dx.doi.org/10.1063/1. 
2711806. 

[25] T. Tarn, D. Dong, and J. Lam, lET Control Theory 
& Applications 3, 161 (2009), ISSN 1751-8652, URL 
http://dx.doi.org/10.1049/iet-cta:20080047. 

[26] K. Jacobs, New Journal of Physics 12, 043005 (2010), 
ISSN 1367-2630, URL http://dx.doi.org/10.1088/ 
1367-2630/12/4/043005. 

[27] M. K. Pedersen, J. J. W. H. Sprensen, M. C. Tichy, 
and J. F. Sherson, New J. Phys. 16, 113038 (2014), 
ISSN 1367-2630, URL http://dx.doi.org/10.1088/ 
1367-2630/16/11/113038. 

[28] S. Ashhab and F. Nori, Phys. Rev. A 82 (2010), 
ISSN 1094-1622, URL http://dx.doi.org/10.1103/ 
PhysRevA.82.062103. 

[29] S. Fu, G. Shi, A. Proutiere, and M. R. James, Phys. Rev. 
A 90 (2014), ISSN 1094-1622, URL http: //dx.doi.org/ 
10.1103/PhysRevA.90.062328. 

[30] Y. Wang, R. Wu, X. Chen, Y. Ge, J. Shi, H. Ra¬ 
bitz, and F. Shuang, J Math Ghem 49, 507 (2010), 
ISSN 1572-8897, URL http://dx.doi.org/10.1007/ 
S10910-010-9755-6. 


